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The three Comments by Barclay and Maxwell [Q, Duncan and Willey and Samuel 

and Steinfelds all correctly point out that my estimate [|[ for r^(l), the nth order 

-\ - coefficient in the perturbative expansion of the normalized total e'^e~ cross-section, does 

^ \ not exhibit the correct flavor or color dependence of the exact calculation Q . This naturally 

casts serious doubts upon the validity of the result. I would like to suggest, however, that, 

^ \ in spite of this, my original estimate remains valid as an asymptotic formula for sufficiently 

large n and that the leading corrections of 0(l/n) have a strong flavor and color dependence. 

Below I give some arguments as to why this might be expected to be so. 

T 1 

^ . It turns out that, if the region of validity of the asymptotic formula is n ^ no{nf, Nc), 

O I then typically no is expected to be relatively large, i.e. > 0(2 — 3). I propose, however, 

^ ! that there is a valley in {nf,Nc) parameter space whose bottom is the approximate line 

CN ■ ^/ ~ 2A^c — 1 and where no^O(l). This passes through, or near, the point (nj = 5, A^c = 3) 

so, in this sense, the close agreement of my result with the exact calculation for the physical 

p. 



m 



case of interest is, indeed, fortuitous. In the general case of arbitrary Uf and A''^ one would 
^ I therefore expect to have to go to much larger values of n (>6, say) to obtain reasonable 



X 



agreement. 

The original point of my paper was to derive an asymptotic estimate for r^ that gave 
the correct order of magnitude: speciflcally to answer the question "is r3 ~ 5 or 50?" This 
was in response to the confusion resulting from the incorrect "exact" calculation which 
originally gave r^ ~ 70. Even though corrections to my estimate were discussed in my 
paper, no serious attempt to evaluate them was made. Since the position of the leading 
saddle point is at ki ~ bi[{n — 1) + 6'], it was natural to retain the combination (n + b') in 
the final result. However as I pointed out, this, of course, does not incorporate all 0(l/n) 
corrections. There are several other sources of 1/n contributions such as corrections in 
going from the dn to the r„ [see my eq. (12)] and from approximating ImD by its leading 
term [see eq. (11)]. The former leads to contributions to r^ like {Sn'^/3){b2/bi)r2 and 
(47r^) (63/361) whereas the latter gives an overall modifying factor [H-l/3{6'4-r2/27r^6i}]. 
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To these must be added corrections to the Gaussian approximation of the saddle-point 
integration. Although a careful systematic examination of all the corrections has not yet 
been carried out (it is presently underway) it is clear that they are, in general, large and 



have a strong dependence on n/ and A''^. For example, in MS with A''^ = 3 and Uf = 5, the 
factor r2/2n'^bi ~ 1.5, whereas with Nc = 5 and n/ = 1, it is almost 3. Furthermore the 
contribution {87i'^/3){b2/bi)r2 which is 1.18 for Nc = 3,nf = 5, is 14 when Nc = 5, Nf = 1! 
Similarly, the term (47r^)^63/36i is only 0.5 for A'c = 3, n/ = 5 but is 7 at A'c = 5, n/ = 1. 
A preliminary evaluation of this set of corrections indicates that, in general, they can be 
uncontrollably large in some cases. However, when Nc = 3 and Uf = 5 they are relatively 
small so that use of the asymptotic estimate for r^ can be justified. On the other hand, 
the remarkable closeness to the exact result is clearly accidental. 

Notice, incidentally, that scheme dependence enters via these non-leading contribu- 
tions. The fact that the leading term is scheme invariant is not an argument against its 
validity. On the contrary, one can argue on very general grounds that the leading large n 
behavior of r„ should, in fact, be scheme invariant. The point is that this behavior deter- 
mines the nature of the divergence of the perturbation series which is itself a refiection of 
the singularity structure in (7^ near (7^ = 0. However, the analytic structure in (7^ can be 
determined via the renormalization group (RG) since this requires that q'^ and (7^ always 
occur in the combination q'^ e^ ^^' [withK (g) = J -gf-r] and the analytic properties in q^ are 
known [see my eq. (9)]. Using the perturbative expansion for P{g) around g^ = gives 
K{g) ^ 1/61(7^ + b' In (7^ + 0{g'^) when (7^ ^ 0. The neglected terms are analytic at g^ = 0. 
The non-analytic structure at g^ = is therefore completely determined by 61 and 62 both 
of which are scheme invariant. This therefore shows (i) that the leading large n behavior 
of the Vn can, in principle, be determined from the RG and q^ analyticity and (ii) that the 
result will depend only on hi and 62 and therefore be scheme-invariant. 

As a corollary, this also demonstrates the importance of 62 since its presence dramat- 
ically changes the analytic structure in (7^. From the fact that there are discontinuities 
only along the positive real axis in q^ one can deduce [see my eqn. (9)] that the g^ (or 
k = 1/(7^) singularities occur when k/bi + b'ln{k+ f-) + - ■ ■ = lnz±2'n:Ni with A an integer 
and z running from zero to infinity. When 62 = this implies that the /c-plane [where the 
integration is to be performed] separates into an infinite number of disconnected sectors 
parallel to the real axis each separated from the next by 2nibi. The appropriate region 
of integration therefore reduces to —00 <Rek< + 00 and < Inik < 2nbi, the boundary 
being the integration contour. This invalidates a derivation of the null result for dn claimed 



for this case (62 = 0) in [|T|]. However, it does emphasize a point that was suppressed in 
my paper, namely that great care must be taken in defining the contour and domain of 
integration as determined by the RG before interchanging the z and k integrals. Indeed, 
when y = 0, this results in no n! behavior from the saddle-point integration. However, if 
y is now included, additional singularities are present, for, when z ^ 0, one can now have 
k ~ —6'. The structure in the A;-plane is now quite different and leads, via a saddle-point 
integration to the n! growth of r^ and the result quoted in my paper. The presence of 62 
is crucial; even though it plays a minor role in the final expression, the result is simply not 
deriveable without it. Thus, conclusions based on its omission, such as in ref. 1 and in the 
work of Brown and Yaffe are not directly relevant. 

The problem of interchanging integrations, which requires care in defining the domain 
and contour of integration, is nicely illustrated by the example of eq. (2) in ref. 1. To 
avoid the "zero times infinity" problem that can occur from a cavalier interchange, the 
(7^-contour in this case should first be wrapped around the cut associated with (— (7^)^~* 
to obtain 

d(s) = sin7r. /■ ^(/)-i-I)(/) (1) 



with L lying above the cut on the positive real axis. This is the appropriate coefficient 
generating function for this case. Now, when the representation eq. (1) or ref. 1 is inserted, 
no problem arises upon interchange of integrations; thus 



dis) = '-^^^ dzfiz) dg'{gY'-'e-/^' 

^ Jo Jo 

f^ dz 
= r{s){smns / —f{z)z-'}. 
Jo TT 



(2) 



The quantity in curly brackets is just the sth coefficient in the perturbative expansion of 
f{z) as it must be. 

Many of these points will be expanded upon in greater detail in forthcoming papers. 



References 

[1] D. T. Barclay and C. J. Maxwell, Phys. Rev. Letters to be published. 

[2] A. Duncan and R. Willey ibid. 

[3] Mark A. Samuel and Eric Steinfelds, ibid. 

[4] G. B. West, ibid 67, 1388 (1991). The notation and equation numbers in this note all 
follow those in this reference. 

[5] L. R. Surguladze and M. A. Samuel, Phys. Rev. Lett. 66, 560 (1991); S. G. Gorishny, 
A. L. Kataev and S. A. Larin, Phys. Lett B259, 144 (1991). 

[6] S. Gorishny, A. L. Kataev and S. A. Larin, Phys. Lett. B212, 238 (1988). 

[7] L. S. Brown and L. Yaffe, Phys. Rev. 45 R398 (1992). These authors set 62 = and 
argue that since they cannot derive my result by their method, then analyticity and 
the RG are not sufficient to determine the large n behavior of r^. However, they do 
not actually use the full analyticity properties of D in the full complex q^ plane, only 
the fact that its absorptive part along the positive real axis is determined by R. 

[8] G. B. West, in "Radiative Corrections" edited by N. Dombey and F. Bondjema 
(Plenum, N.Y., 1990) p. 487. 



